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ABSTRACT 

High-dispersion spectroscopic observations are used to refine orbital elements 
for the symbiotic binary CH Cyg. The current radial velocities, added to a 
previously published 13 year time series of infrared velocities for the M giant in 
the CH Cyg symbiotic system, more than double the length of the time series 
to 29 years. The two previously identified velocity periods are confirmed. The 
long period, revised to 15.6 ±0.1 yr, is shown to result from a binary orbit with 
a 0.7 Mq white dwarf and 2 Mq M giant. Mass transfer to the white dwarf is 
responsible for the symbiotic classification. CH Cyg is the longest period S-type 
symbiotic known. Similarities with the longer period D-type systems are noted. 
The 2.1 year period is shown to be on Wood's sequence D, which contains stars 
identified as having long secondary periods (LSP). The cause of the LSP variation 
in CH Cyg and other stars is unknown. From our review of possible causes, we 
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identify g-mode non-radial pulsation as the leading mechanism for LSP variation 
in CH Cyg. If g-mode pulsation is the cause of the LSPs a radiative region is 
required near the photosphere of pulsating AGB stars. 

Subject headings: infrared: stars — binaries: symbiotic — stars:individual (CH 
Cyg) — stars:late-type — stars: variables: other 



INTRODUCTION 



Symbiotic stars are mass transfer binaries, containing a mass losing red giant and a mass 
accre ting companion th at is usually a white dwarf. One system is known with a neutron 
star (IHinkle et al.ll2006l ). and systems containing main sequence stars are possible but have 
proved elusive to identify. Mass transfer causes the symbiotics to exhibit hi ghly complex light 
variat ions and spectra. From their characteristics at infrared wavelengths IWebster fc Allen 
( I1975I ) separate the symbiotics into two subclasses, D for dusty-type and S for stellar-type 
systems. CH Cyg is the brightest of the S-type symbiotic stars. 

The current series of papers is concerned with infrared spectroscopy of the late-type star 
in each symbiotic system. Because velocities of this component are usually well behaved, we 
employ these spectra to derive single-lined spectroscopic orbits. The orbital elements provide 
insight into the period distribution, mass ratios, and orbital dynamics of these interacting 
systems. These results are ultima tely of inter est in evaluating the potential for catastrophic 
demise of the stellar components (jibeni 120031 ). 



CH Cyg was th e topic of our first paper on the infrared spectroscopy of symbiotic stars 
( iHinkle et al.l Il993l . hereafter Paper I). As a result of its brightness and placement in the 
northern sky, CH Cyg is easily observed. Following the publication of Paper I, we continued 
to monitor CH Cyg with various spectrographs and now have data spanning a period of nearly 
three decades. In the intervening 16 years since the publication of Paper I our knowledge 
of M giant variables as well as the CH Cyg system itself has greatly in creased. We find, in 
agreement with many others who have worked on symbiotics (see e.g., ISkopal et al.lll996al : 
Crocker et al.l l2002l : ISokoloski fc KenyonI l2003bl ) , that CH Cyg is the most complex of the 
symbiotics that we have observed. 

Based on velocity variations of the M HI star in the CH Cyg system, we concluded in 
Paper I that the CH Cyg system is triple. Two separate stable velocity variations, a long 
period variation of ~14.5 years and a short period variation of ~2.1 years, were found and 
are confirmed and refined in the present paper. In the discussion section of Paper I we 
proposed a model for this velocity behavior that placed the symbiotic binary in the 2.1 year 
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orbit. Our model, based on the short period mass function, had three underpinnings: 



(1) The absence of any other known S-type symbiotic with a confirmed orbital period 
l onger than five ye ars and a mean period distribution for S-type symbiotic orbits of ~2 years 
jPekel et aliboOTh . 



(2) Stellar pulsation theory showi ng that the two year ve locity variation is too long to 
be the M giant fundamental pulsation (IHughes fc Woodlll990l ). 



(3) Weak evidence at that time that the system had a high inclination. 

Because this triple-star model has been controversial and new information has accu- 
mulated, we decided to reexamine this conclusion. In the present paper we review all the 
evidence and show convincingly that the triple-star model with the white dwarf as a member 
of the short period system is incorrect. We now conclude instead that the system is a binary 
with the mass accreting white dwarf in the long period orbit. 

Nonetheless, points (1) and (2) remain valid and are essential to an understanding of 
CH Cyg. While the CH Cyg symbiotic system has unique properties, we believe the most 
interesting aspect of CH Cyg is the nature of the 2.1 year variation. We show that CH Cyg 
qualifies as the most intensively investigated variable star with long secondary period (LSP) 
light and velo city var i atioii s. The existence and uncertain nature of the LSP variations is 
highlighted bv lWoodI kQQl\ }. 



OBSERVATIONS AND REDUCTIONS 



In 1992 we obtained two spectroscopic observations with the Pourier transform spec- 
trometer ( PTS) at the co ude focus of the Kitt Peak National Observatory (KPNO) 4 m 
telescope ( iHall et al.lll978l ). These observations are a continuation of the series of 71 PTS 
spectra discussed in detail in Paper I. The two observations were obtained in the 2 fim region 
and have a resolving power of ~60,000. Shortly after the two new spectra were obtained, 
the PTS was decommissioned as a result of budget cuts at KPNO. 

Prom 1995 through 2000 we collected 25 observations with the coude feed telescope and 
spectrograph system at KPNO. The detector was an infrared camera, NICMASS, developed 
at the University of Massachusetts. We obtained a 2 pixel resolving power of 44,000 at a 
w ayelength of 1 . 623 At m. A more extens i ye des cription of the experimental setup can be found 



m 



Jovce et all fll998f ) and IPekel et al.l (120001 ). The NICMASS equipment was sent to Mt. 



Stromlo Observatory to continue the symbiotic orbit program in the southern hemisphere. 
That equipment along with the 1.85 m telescope was destroyed in a devastating bush fire 
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that burned over Mt. Stromlo in January 2003. 



In 2000 three spectra were also acquired with the Phoenix cryogenic echelle spectrograph 
at the f/15 Cassegrain focus of either the KPNO 2.1 or 4 m telescopes. A complete description 
of the spectrograph can be found in lHinkle et al.l (119981 ). The observations were centered at 
either 1.563 or 2.226 /im and have resolvin g powers of either 50,000 or 70,000. An expanded 
discussion of the experimental setup is in iFekel et al.l (120001 ). Phoenix was deployed to the 
southern hemisphere in 2000. 

Finally, five spectrograms were obtained in 2007 and 2008 with the KPNO coude feed 
telescope, coude spectrograph, and a CCD with enhanced red sensitivity, designated LBIA. 
This 1980 X 800 pixel CCD was manufactured by Lawrence Berkeley National Laboratory 
and is 300 /xm thick. Although this thickness results in increased pixel contamination by 
cosmic-ray and background events, the chip was used because of its high quantum efficiency 
at far-red wavelengths. Our spectrograms, centered near 1.005 /im, have a wavelength range 
of 420 A and a resolving power of ~21,500. 

Reduction and radial velocity measurement of the FTS spectra are discussed in detail in 
Paper I. F or the NICM ASS and Phoenix data standard observing and reduction techniques 
were used ( |Joycelll992l ). Wavelength calibration at the infrared wavelengths of 1.563, 1.623, 
and 2.226 /xm posed a challenge, because the spectral coverage was far too small to include 
a sufficient number of ThAr emission lines for a dispersion solution. Thus, our approach was 
to utilize the absorption lines of a K III star. Several sets of lines were tried, including CO, 
Fe I, and Ti I. These groups all gave consistent results. For the spectrograms acquired with 
the LBIA CCD at 1.005 /im, we were able to use ThAr spectra for wavelength calibration. 
Telluric lines are present in the 1.005 and 2.226 /im wavelength regions. These lines were 
removed from our observations by ratioing the spectra to a hot star spectrum observed on 
the same night. 

For all of our spectra except those obtained with the FTS, the absorption line ra- 
dial velocities of C H Cyg were measured with the IRAF cross-correlation program FXCOR 
(iFitzpatrickl Il993l ) . Those velocities were determined relative to the M-giant International 
Astronomical Union velocity sta ndards 6 Oph or a Cet, which have radial velocities of —19.1 
and —25.3 km s~\ respectively (IScarfe et al.lll990l ). The standards were observed multiple 
times during the course of each night. For the FTS spectra, which cover a wavelength range 
that is ~100 times greater than the other spectra, velocities were referenced to telluric lines. 
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3. ORBITAL SOLUTION 

Our 106 KPNO observations, which span nearly 30 years, are hsted in Table 1, and the 
velocities are plotted in Figure 1 versus hehocentric Julian date. That plot cl early sho A vs two 



periodicities. Thus, we have employed the general least squares program of iDanielsl (Il966l ) 
to obtain a simultaneous orbital solution for the short- and long-period velocity variations of 
CH Cyg. All KPNO velocities were given unit weight after a comparison of the FTS velocities 
with those from the Phoenix and coude feed spectrographs showed that the RMS values of 
the velocity data sets were nearly identical. We then computed several different orbital 
solutions. The first used only the KP NO velocities, while, s imilar to a solution of Paper I, 



the second included the velo cities of iDeutsch et al.l (119741 ). lYamashita fc Maeharal (119791 ) 



and lYoo &: Yamashital (Il99ll ) plus our KPNO velocities and is designated the all-velocity 
solution. Weights adopted for the non-KPNO velocities are the same as those assigned in 
Paper I. As was done in that paper, we have also computed orbital solutions with the short 
period eccentricity fixed at 0.0. However, given the systematic residuals that result from 
adopting such a sinusoidal fit and our discussion and conclusion about whether this short 
period velocity variation results from orbital motion or pulsation, we do not discuss the e = 
0.0 solutions any further. 

A comparison of the two remaining solutions shows that the long period of the KPNO- 
only solution is nearly 40 days longer than that of the all- velocity solution, but that difference 
is only 0.7% and is less than the period uncertainty. Uncertainties of the orbital parameters 
in the all-velocity solution are somewhat smaller than those in the KPNO-only solution, as 
would be expected because of the longer baseline of data, but the parameter values of both 
solutions are within those uncertainties. Thus, in Table 2 we have chosen to adopt the more 
homogeneous solution that includes only our KPNO velocities. 

For the individual observations Table 1 lists the heliocentric Julian date, the observed 
total velocity, and the observed minus calculated velocity residual (O — C) to the combined 
orbit. Also computed and listed in the table are the long period orbital phase, the long period 
velocity, which is equal to the total velocity minus the computed short period velocity, the 
short period orbital phase, and finally, the short period velocity, which is equal to the total 
velocity minus the computed long period velocity. Figure 2 presents the computed velocity 
curve of the long period orbit compared with the KPNO radial velocities, where zero phase 
is a time of periastron. Each plotted velocity consists of the total observed velocity minus 
its calculated short period velocity. Figure 3 shows the computed velocity curve of the 
short period orbit compared with the KPNO radial velocities, where zero phase is a time of 
periastron. Each plotted velocity consists of the total observed velocity minus its calculated 
long period velocity. 
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The radial velocity residuals, listed in column 3 of Table 1 ("O — C"), were analyzed 
for possible additional periods. The program PeriodoGRAM (PGRAM) was employed to 
fit sinusoids to the phased data for periods between 50 and 500 days in steps of 0.05 days. 
This period finding program works well even for non-sinusoidal curves with eccentricities as 
large as ~0.4. Various subsets of data were searched in an attempt to remove false periods. 
Two possible periods of ~62 and ~170 days appear in the velocity residuals. However, the 
velocity semiamplitudes are very small, at maximum 0.4 km s~^. Of particular interest, 
we find no period above the noise in the 90-110 day period range, which corresponds to a 
photometric period of CH Cyg. 



Yamashita fc Maeharal (Il979l ) first proposed an orbit for the M giant. Despite the use of 
very heterogeneous data with large uncertainties, they determined a probable period of 5750 
days, very close to our value of 5689 days. Although their semiamphtude and eccentricity are 
rather different from ours, other elements are similar. We also compare our orbital solution 
to those in Paper I. In that earlier analysis the long period for the all-velocity solution was 
5298 days, while that for the KPNO velocities alone was 5483 days. Thus, the improved 
period of our adopted solution (Table 2) is nearly 400 days longer than the all-velocity 
solution in Paper I and about 200 days longer than the KPNO-only solution given in that 
paper. Comparing other elements of our adopted solution and the KPNO-only solution of 
Paper I, the orbital elements are similar, but the uncertainties of the long period elements 
have been significantly reduced. 

The ephemeris for conjunctions with the M giant in front of the white dwarf, which 
corresponds to times of mid-eclipse, is 

T^onjiHJD) = 2, 446, 353(±192) + 5689(±47)^, (1) 

where E is an integer number of cycles. The above ephemeris for the eclipse of the white 
dwarf by the M giant predicts mid-eclipse dates of HJD 2,440,664 for 1 cycle earlier and 
HJD 2,452,042 for 1 cycle later. The uncertainty of these predictions is about 240 days. 



4. BASIC SYSTEM PARAMETERS 

Paper I presented a broad introduction to CH Cyg, but much material has been pub- 
lished subsequent to that time. To understand the nature of the CH Cyg system, we first 
undertake an extensive review of the literature to determine basic parameters: variability, 
distance, luminosity, radius, rotation, mass, mass loss, and orbital inclination. This results 
in a revised picture of the system. 
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4.1. Variability 



The photometric va r iability of CH Cyg was reviewed by iMuciek fc Mikolajewskil (119891 ) 
and iMikolajewski et al.l (Il990al ). The mass accreting component in the CH Cyg system 
contributes to the flux especially at blue wavelengths. During periods of strong activity, the 
M giant variability can be cloaked even in the optical with light generated by the accretion 
process. However, during quiescent periods, the M giant can be detected into the blue. 

CH Cyg had an extended period of quiescence from 1885 - 1963. During this interval CH 
Cyg was known only as a late-t ype semiregular variable star. A review of historic material 
on the variability of CH Cyg by lMuciek &: Mikolajewskil ( 1l989l ). obtained during the period 
of quiescence, flnds that the dominant period in the M giant is ~100 days. During pe- 



riods of symbiotic activit y the 100 day period can be hard to detect flMikolajewski et al. 



1992: Munari et al 



199 



The visual a mplitude of the 100 day period is small, ~0.1 



magni tude (jMuciek fc Mikolaiewskil 119891 ). The exact period depends on the dates sam- 
pled. [Mikolaiewski et al. (1990a) gave the dominant period as 94 an d 99±2 days while 
Mikolaiewski et all Jl992h found 102±3 days. IMikolajewski et al.l (119921 ) concluded that the 
~100 day period is sy stematically inc r easing at a rate of ~5 days per century and suggested, 
following the work of IWood fc Zarrol (Il981l). t h at the CH Cyg M giant is curr e ntly u nder- 
going a helium-shell flash. IMikolajewski et al.l (Il990al ) and IMikolajewski et al.l (119921 ) both 
noted that the ~100 day period implies that the CH Cyg giant is a flrst overtone pulsator. 

In addition to the 10 day period, a period of ~770 days is also frequently reported. 
Mikolajewski et al.l (Il992l ) provided a detailed analysis of the 770 day photometric period 
in CH Cyg. The photometric properties of the 770 day variation show increasing amplitude 
toward shorter wavelengths but not in a way consistent with dust extinction. In addition, 
the energy distribution in the IR suggests changing spectral type (and tempe rature) during 



the lig ht cycle, but there are no s pectral type c hange s in the blue. While iMunari et al. 



( 1l996l ) failed to detect this period, ISkopal et al.l (120071 ) noted the presence of a ~750 day 
period with amplitu de of nearly 1 mag i n V photometry. The ~750 day period is apparent 
in data illustrated in lSkopal et al.l (120071 ). The symbiotic was quiescent at the time of these 
observations. There also is a suggestion of the ~100 day period in this data, but the spacing 
of the data is too coarse for the period to be adequately sampled. 



The simultaneous presence of both periods in CH Cyg was reported by lPayne-Gaposhkin 



( 119541 ). She listed CH Cyg as one of a family of pulsators with a LSP and estimated a 
typical long to short period ratio of ~9 for the s e obj ects. A longer list of LSP semiregular 
variables appears i n iHoukI (Il963l). IWood et al.l (120041 ) asserted that the LSP stars listed in 
Payne-GaposhkinI (119541) an d iHoukI (119631 ) are those with the largest amplitude secondary 
periods. The iHoukI (119631 ) catalog contains only ~1.5% of the then known long period 
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variables. More rece ntly, the fraction o f loc al long period vari ables with long secondary 



periods is reported by IPercy et al.l (12004 ) and IWood et al.l (120041 ) as ~25-30% 



In paper I we noted the close agreement between the 770 day photometric period and the 
750 day spectroscopic period and assumed bot h originate fro r n the s ame physical mechanisn i. 
We make the sa r ne ass umption in this paper. iHinkle et al.l (120021 ). lOlivier fc WoodI (120031 ). 
and lWood et al.l (l2004l ) detected LSP variations spectroscopically in a number of semiregular 
variables. The semiregular 100 day pulsation was not detected spectrosco pically. Scaling 



from typical ratios of velocity to photometric amplitude for SRb variables (iLebzelter et al. 
2000), the CH Cyg 0.1 magnitude visual amplitude implies a ~0.5 km s~^ velocity amplitude. 
This is near the detection limit for our observations (§3). 



4-1.1. Eclipses 



When the accretion disk in the CH Cyg system is active, light from the disk dominates 
the blue a nd violet por t ions o f the spectrum. At such times the optical light curve is highly 
complex. lEyres et al.l (120021 ) present a summary of recent optical photometry. Because 
CH Cyg has a complex light curve, evidence for eclipses depends critically on collaborating 
evidence rather than just sudden decreases in the U band light. The first claim of an eclipse 
of the hot component by the red giant, occurring from 1985 May to October, was reported 
by lMikolajewski et al.l (119871 ). This was largely based on timing, matching a U band event to 
one more than 15 years earlier. However, these decreases in light show considerable structure 
because of activity in the system. The dates of mid-eclipse are JD 24405 85 and JD 2446270 
with th e eclipse FWHM ~200 days. The ingress is brief, perhaps 10 days. iMikolajewski et al. 
Jl990ah reported an ephemeris of JD^i„ = 2,446,275(±75) + 5700(±75)E. 



Subsequent to the claim of IMikolajewski et al.l (119871 ) considerable evide nce has been 



publish ed that supports the claim that the 1985 U band decrease is an eclipse. ISkopal et al. 



(Il996al ) summarized a list of five changes in the CH Cyg system that were observed during the 
1985 eclipse. These are (1) the disappearance of rapid optical flickering, (2) the dominance of 
the M giant continuum in optical spectra at mid-eclipse, (3) the transition of double peaked 
Balmer emission lines to weak broad single peaked emission lines at mid-eclipse, (4) the 
dominance of the red side of the double peaked emission in Ha and H/3 before mid-eclipse 
and the dominance of the blue emission following mid-eclipse, and (5) the prominence of the 
nebular lines [Ne III] 3869 A and [O HI] 5007 A. Items (1) and (2) show that the accretor 
was eclipsed. Items (3) and (4) identify the eclipsed object as a rotating disk. Item (5) 
shows that the blue - ultraviolet flux level was depressed during the eclipse so that nebular 
lines became visible. 
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Apparent confirming proof of eclipses in. CH Cyg comes from tfie predict i on and tfien 



observ ation of the eclipse of 1999 (lEyres et al.ll2002l : ISokoloski &: Kenyonll2003bl ) . lEyres et al. 



( 120021 ) presented UBV photometry for the 1999 eclipse, revealing that the eclipse is not 
visible at V, but is clearly detectable in B and much more pronounced in U. This is the 
expected enhancement of an eclipse of a hot object by a cool star. In addition, they find that 
flickering is again absent during the 1999 eclipse, just as in 1985. The strengt hs of various 



accretion related spectral lines decreased during the eclipse. From Figure 5 of ISkopal et al. 
(120071 ) the mid-date of the eclipse in 1999 is estimated as JD 2,451,425. 



A major complication to this recent eclipse identification is t hat the 1999 eclipse date, JD 
2,451,425, is 550 days earlier than predicted by the ephemeris of iMikolajewski et al.l (Il990al ) 
and 617 days earlier than our ephemeris. These differences correspond to an orbital phase 
shift of ~0.1. In eclipses of symbiotic stars the observed eclipse is not that of the white dwarf, 
which is unobservable at visual wavelengths, but rather of a UV bri ght, hot spot. Three 
dimensional simulations of the symbiotic-recurrent nova RS Oph by IWalder et al.l (120081 ) 
predict an accretion disk that has a diameter 0.1 times the major axis of the binary system. 
This accretion disk is surrounded by a much larger Archimedean spiral, so ther e is a lar g e area 
over which the hot spot can occur. From observations of several symbiotics ISkopall (119981 ) 
concluded that during quiescent phases of the systems, the times of eclipses occur prior to 
the time of inferior conjunction of the giant. When the systems are acti ve eclipses coincide 



with co ni unctions predicted by spectroscopic orbits. For CH Cyg the iMikolajewski et al. 



(jl990al ) ephemeris from active periods agrees with the times of inferior conjunction given in 
Equation 1. The 1999 eclipse took place during the period of declining activ ity preceding 
the current extended period of quiescence (lEyres et al.ll2002l : ISkopal et al.ll2007l ). We suggest 
that the accretion disk was undergoing changes at the time of the 1999 eclipse. 

The details of the light curve shape and hydrogen line profiles in these eclipses add 
further support to the eclipse interpretation. In all cases the U light is brighter before the 
eclipse than after. In cataclysmic systems this is a well-known eclipse hallmark. The pre- 
eclipse brightening results from the accreting side of the disk facing the mass donating star 
and rot ating in the direction of orbital motion . For Ha and H/3 line profiles, observed through 
eclipse ( lEyres et al.l 120021 : iFernie et al.lll986l ). the negative side of the profile is eclipsed first 
and reappears first. Thus the accretion disk is rotating with the side moving away from the 
red giant on the leading side of the orbit. 

The eclipses discussed above are associated with the 15 . 6 yr o r bit. Yet anot h er con i- 
plication in the discussion of CH Cyg eclipses is that ISkopall (Il995l ). ISkopal et al.l (Il996al ). 
and su b seque nt papers by Skopal claim that there are also eclipses in the 2.1 yr orbit. 
Skopall (119951 ) reported the presence of [/-band decreases spaced at ~2.1 year intervals on 
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the light curve. Supporting evidence is seen in the disappear ance of flickering d u ring at 
leas t some of these events. Cou nter to the 2.1 yr eclipse claim iMikolajewski et al.l (Il990bl ) 
and ISokoloski &: Kenyoru (j2003bl ) found a general co rrelation of the streng t h of th e flickering 



with the luminosity of the blue-violet continuum . ISokoloski fc KenyonI (l2003al ) suggested 
that the events observed by ISkopal et al.l (Il996al ) are not eclipses but result from collapses 
of the inner accretion disk. 



4.1.2. Flickering and Jets 



Flickering is an indicator of a compact mass-accreting object in the system. IMikolajewski et al. 



(1l990bl ) found that the rapid optical variations (flickering) of CH Cyg have a pe riod of ~500 s . 
The re is a range in the pe riods reported for optical flickering in CH Cyg, e.g.. iHoardI (119931 ) 
and Rodgers et al. (Il997) suggested a pe riod in the range ~ 2200 - 3000 s (for a detailed 



discussion see ISokoloski fc KenyonI l2003bl ). With any of these time scales the flickering must 
originate at a hot component in the CH Cyg system since the dynamical time scale for 
the giant, trfj^„~(i?'^/2G'M)^/^, is on the order of a month. Flicke ring is also seen in X-ray 
observations, conflrming the white dwarf nature of the secondary (lEzuka et al.lll998l ). 



CH Cyg is known to have both radio (ITaylor et al.lll986l ) and optical (ISolflll987l ) jets 
which are approximately in the plane of the sky. This provides a flrst limit on the inclination 



of the orbit, although the information from the eclipses is far more accurate. ICrocker et al. 



( I2OO2I ) found from a time series of observations that the jet is precessing with a period ~18 ± 
0.5 yr. While the precession mechanism is uncertain, the period of the precession is similar 
to the 15.6 yr orbit, suggesting that the white dwarf is in that orbit. 



Sokoloski fc KenyonI (l2003al ) discussed the relation between the inner disk and the devel- 



opment of jets. They argued that it is impossible to understand the jet and optical activity 
of CH C yg without an accretion disk around a white dwarf companion. ISokoloski &: Kenyon 
( l2003al ) pointed out a strong correlation between optical flickering in CH Cyg and changes 
in flux. The decrease or disappearance of flickering corresponds to a decrease in the blue 
flux. They found that this is related to disruption of the inner disk and follows episodes of 
jet production. 
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4.2. Distance, Luminosity, and Radius 



4.2.1. Giant 



Kenyon fc Fernandez-Castrd (119871 ) undertook a detailed photometric c alibration of 
symbi otic spectral types. From the resulting type of M6.5 (± 0.3) 111 for CH Cyg. lMiirset et al, 
( 19911 ) used a spectroscopic parallax to derive a distance of 240 P*^- "^^^^ excellent 



agreement with the 268±66 pc distance based on Hipparcos observation s (jViotti et al.lll997l) . 
which has been used extensively in the literature since 1997. Recently, Ivan LeeuwenI (120071 ) 
reanalyzed the Hipparcos parallaxe s resulting in a dis tance for CH Cyg that is nearly 10 



percent closer, 244 pc, than the IViotti et al.l (119971 ) distance. The uncertainties for the 



original and revised Hipparcos distances have large overlap. 



Biller et al.l (120061 ) provided a fit to optical and infrared photometry of CH Cyg, using 



a Kurucz model spectrum for an M6 HI star and a dust model. They found a l uminosity 
of 6900 Lq. This is in good agreement with the value given by ISkopall (119971 ) from an 
energy distribution encompassing the UV through the far-lR. He determined a bolometric 
luminosity of ~8000 ± 4000 Lq for a distance of 268 ± 66 pc. This luminosity corresponds 
to a radius of ~310 R©. 



Dyck et al.l (119981 ) measured a 10.4 mas diam eter for CH Cyg at 2.2 /im. Assuming the 



244 pc distance, the stellar radius is then 2 73 R^. ISchi 



65 Rq , based on near-IR photometry. The iDyck et al 



d et all (Il999f ) find a radius of 280 ± 



(119981 ) radius can als o be calibrated 



to an effective temperature of 3084 ± 130 K. This compares very well with the lRichichi et al. 
(Il999h calibration for spectral type M 7 of 3150 ± 95 K. 



CH Cyg has a variable K mag ranging from ~— 0.9 to —0.3 (iMunari et al.lll996f). We 
adopt the mean as a typical magnitude, K = —0.6 and J — K = 1.6. iTaranova &: Shenavrin 
(I2OO7I ) reported K band fading, with K as faint as +0.35 and J — K = 2.06, but this appears 
to be related to dust formation episodes. With the 244 pc distance, K = —0.6 corresponds to 
an absolute K magnitude of —7. 5 =fc 0.4. The K band b olometric correction can be computed 
from the J — K color, following iBessel fc WoodI (119841 ). to yield a bolometric magnitude of 



—4.3 and a luminosity of 4200 ^^^^^ Lq. Assuming a temperature of 3100 K, the effective 
temperature - radius - luminosity relation gives a radius of ~224 R©. 

The values of the M giant's basic parameters, derived from the various analyses discussed 
above, are in reasonable agreement. For the following discussion we adopt values of L ~5000 
Lq, T ~3100 K, and R ~280 Rq. The corresponding bolometric magnitude is —4.5. We 
note that typi cal radii for other giant stars in symbiotic binaries of earlier spectral types are 
100 - 200 Rq jFekel et al.ll2003l ). 
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4.2.2. White Dwarf 



While the white dwarf in the CH Cyg system has not been observe d directly, some 



prope rties can be inferred from the accretion process. As summarized by iKarovska et al. 



( 120071 ). X-rays from CH Cyg have been detected by a number of groups. The X-rays originate 
from the accretion disk - white dwarf bo undary layer a s wel l as from a shock, where the 
jet interacts with circumbinary material. lEzuka et al.l (119981 ) uses X-ray spectroscopy to 
determine the temperature and bolometric luminosity of the white dwarf. The luminosity 
derived, ~10^^ ergs s^^ (~ 0.25 Lq), can be combined with a mass-radius relation to set a 
lower limit on the mass of the white dwarf of 0.44 Mq. 



4.3. Rotation and Line Widths 



Schmutz et al.l (Il994j ) , iMiirset et al.l (120001 ) , and IZamanov et al.l (120071 ) argue that in 
most S-type symbiotics the giant star is synchronously rotating. In the case of synchronous 
rotation the projec ted rotational velo city of the late-type giant can be used to estimate the 
stellar radius (e.g., iFekel et al.ll2008l ). The v sin i value for the M giant can be estimated 
from the full-width at half-maximum (FWHM) of unblended absorption lines in the near- 
i nfrar ed spectra. Using atomic lines near 2.223 /xm, we apply the analysis technique of iFekel 
(119971 ) to CH Cyg. For CH Cyg an empirical calibration with intrinsic line broadening of 3 
km yields v sin z of 8 ± 1 km s~^. 

Both the 15.6 yr and 2.1 yr orbits of CH Cyg are eccentric. In an eccentric orbit the 
rotational angular velocity of the M giant synchronizes with the orbita l angular velocity 
at peri astron. Th is was called "pseudo-synchronization" by iHutI (Il98ll ). Using equation 
(42) of iHutI (1198 ll ). we calculate a long pseudo-synchronous period of 5219 days and a short 
pseudo-synchronous period of 453 days. 

Assuming sin i = 1, the stellar equatorial r otational velocity is 8 km s^^. Combining 
this value with the pseudo-synchronous periods (IFekel et al.ll2003l ) suggests a stellar radius 
of 825 R0 for the 15.6 yr orbit or 72 Rq for the 2.1 yr orbit. The two computed radii 
are in poor agreement with the ~280 Rq value from the direct measurement of the stellar 
angular diameter. Thus in both the long and short period cases the rotation of CH Cyg 
is not pseudo-synchronous. O ther symbiotic sy s tems with orbital pe riods near that of the 
short period are synchronized (IFekel et al.l 120031 : IZamanov et al.l 120071 ). 



In Paper I we found that for CH Cyg lines of modest strength the FWHM varies by ~ 
30% over the 2.1 yr period. The pseudo-synchronous period differs by nearly a factor of two 
from the spectroscopic period. However, the line width modulation is clearly occurring with 
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the 2.1 yr period. The rotational broadening would be more than twice that observed for a 
period of 2.1 yr and a diameter equal to that of CH Cyg. We conclude that the variation in 
the line widths is not due to a rotational modulation associated with star spots or plages. 



4.4. Masses 



The values for the stellar parameters of CH Cyg can b e used with stellar evolution 
theory to place additional constraints on m ass. ISchmidt et al.l ( 20061) d etermine a solar iron 
abundance for CH Cyg. From Table (2) of I Vassiliadis fc WoodI (119931 ) the solar abundance 
and a bolometric magnitude of —4.5 for CH Cyg yield an initial mass for the M giant between 
~1.5 and 3.5 Mq. Other arguments also suggest a mass near 2 Mq. If CH Cyg is currently 
undergoing a helium-shell flash (§4.1), the luminosity is currently near a peak, which requires 
a mass near t he bottom of the ma ss range. Based on the absolute magnitude and effective 
temperature, ISchmidt et al.l (120061 ) find log g ~0 for the surface gravity of CH Cyg. This 
also is in agreement with a mass of 2 for a radius of 280 R©. 

If we assume th at the white d A varf is not the result of a pathological evolutionary process, 
then the relation of iKalirai et al.l (120081 ) can be used to map the white dwarf initial to final 
mass. The initial mass of the white dwarf, which evolved from the more massive star in the 
binary, had to be larger than the initial mass of the current M giant requiring the mass of 
the white dwarf to be >0.56 M©. 



4.5. Mass Loss 



Gas an d dust mass loss rate s for a number of symbiotics, based on IRAS photometry, 
are given by iKenyon et al.l (119881 ). Revising the distance to 244 pc, the dust mass loss rate 
of CH Cyg is 4 X 10~^ Mq yr~^ and the gas mass loss rate is 1 x 10~^ yr~^. A m ass loss 



rate for dust, based on near-IR p hotometry, is computed by iTaranova fc ShenavrinI (120041 ) 
andlTaranova fc ShenavrinI (120071 ). They obtain ~3 x 10"'' yr~^ in good agreement with 
Kenyon et al.l (119881 ). A dust formation episode in late 2006 resulted in multiple magnitude 
dimming in the J, H, and K bands plus dimming o f a few tenths of a magn i tude a t L and M. 
Assuming a spherical symmetry for the dust shell, iTaranova &: ShenavrinI (120071 ) translated 
the brightness change into a mass-loss rate of 2 x 10~^ M0 yr~^. 



Kenyon et al.l (119881 ) der ived a dust temperature of 4 00 K, considerably less than the 



value of 750-800 K given by ITaranova &: ShenavrinI (120041) . The radius of the d ust scales 
inversely to the dust temperature (IKenyon et al.lll988l ). From lKenyon et al.l (119881 ) equation 
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(5), Rd ~ 108 AU for 400 K dust and 19 AU for 800 K dust. The semi-major axis of the 
15.6 yr orbit is on the order of 8 AU, so the dust at either temperature hes well outside the 
stellar orbits. 



Mid-infrared imaging shows a Gaussian dust shell of ~45 AU FWHM (IBiller et al.ll2006l ). 
This confir ms the circumbinary n ature of the dust and sugges t s that the dust temperature is 

ues a nd in the same 
2007h. Because the 



between the Kenyon et al. ( 1988 ) and Taranova &: Shenavrin ( 2004j) 



temperature range as found in D-type systems (IKotnik-Karuza et al. 



va. 



dust is circumbinary, counter to the assumption of iTaranova fc Shenavrin! (120071 ). extinction 
events are no doubt due to localized condensations (clouds) rather than large scale mass 
ejection events. 



4.6. Inclination 



Fekel et al.l (120081 ) note that the minimum inclination for symbiotic eclipses is given by 

cos{i) < {Rrg + Ryod)/a, (2) 

where a is the semimajor axis of the binary, i is the orbital inclination, and R is the radius 
of the red giant and white dwarf components. From §4.2.1, R = 280 R©. From the orbital 
elements the red giant semimajor axis can be determined, but it is the total semimajor axis 
that is required for equation (2), so a total mass must be assumed. For a first estimate we 
adopt the total mass from stellar evolution (§4.4) of 2.5 Mq. Then from Kepler's third law 
and the 15.6 yr period the semimajor axis is 8.5 AU resulting in a minimum inclination of 
81°. With this estimate for the inclination, sin i, which is critical to converting the observed 
mass function into component masses, is equal to 1.0 to better than 2%. 

The length of the eclipse gives more precise inform ation about the inclination. The 
eclipse duration is ~ 200 days (IMikolajewski et al.l 119871 ). Ingress/egress is short compared 
to the length of the eclipse as expected if the size of the eclipsed UV bright spot is much 
smaller than the diameter of the red giant. The 15.6 yr orbit has a small eccentricity of 
0.122, and so for simplicity we assume the orbit is circular. The eclipse duration is then 
3.5% of the orbital period. For a semimajor axis of 8.5 AU the white dwarf moves 404 Kq 
during an eclipse. This value sets a lower limit on the diameter of the red giant. 

We can compute the inclination by comparing the distance moved during the eclipse 
with the 280 R© radius of the giant. From geometry, the white dwarf traverses a path 194 
Rq off the cent er of th e M gi ant, giving an inclination of 84.0°. This is in agreement with 
an estimate by ISkopall (119951 ) of i > 83°, based on the shape of the 1985 eclipse and the 



assumption of Rrg = 200 Rq. 
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4.7. Summary of Basic Parameters 

From the above sections we can distill a profile of the components of CH Cyg. The M 
giant velocities show the existence of both a short period (2.1 yr) and long period (15.6 yr) 
"orbit" . However, only two stars have been identified directly or indirectly in the CH Cyg 
system, an M giant and a white dwarf. The M giant has L ~5000 Lq (M^o; =—4.5), T ~3100 
K, R ~280 Rq, and M ~2 Mq. The white dwarf has M > 0.6 Mq. The M giant has solar 
abundances, is on the thermal pulsing AGB, and is a radial overtone pulsator with a period 
of ~100 days. The CH Cyg long period system has an orbital period of 15.6 years, eclipses, 
and has an inclination of 84°. These data are summarized in Table 3. 



5. MODELS OF THE CH CYG SYSTEM 

Several models have been proposed to explain the observed properties of the CH Cyg 
system. In the following sections we briefiy review and comment on these models. 



5.1. Mikolajewski et al. 1987 



From the rather poor radial velocities available at that time, iMikolajewski et al.l (119871 ) 
concluded that CH Cyg is a binary with a period ~5700 days and a large orbital eccentricity 
of 0.55. Whi le we have confirmed the p eriod, our radial velocities result in a greatly reduced 
eccentricity. IMikolajewski et al.l (119871 ) estimated that the system consists of an M6 giant 
with a mass of ~3 Mq and a white dwarf companion with a mass ~1 Mq. The white dwarf 
is surrounded by an accretion disk formed from the red giant wind. Eclipses were identified 
in U band light curves, and as a result, the inclination was set to 90°. 



5.2. Hinkle et al. 1993 



Hinkle et al.l (119931 . Paper I) demonstrated that the CH Cyg system has two independent 
sets of velocity variations, a 2.1 year short period and an ~ 15 year long period. They note 
that other well known symbiotic systems have orbits of about 2 years and that CH Cyg can 
exhibit large amounts of activity associated with mass transfer to the secondary. This was 
presumed to imply a small orbital separation. In addition, at this time (early 1990s) the 
evidence for eclipses was weak. To devise a model similar to other active symbiotic systems, 
it was proposed that the components of the symbiotic system reside in the short period orbit. 
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Paper I then used the mass function (iRussell et al.lll955l ) 



/(m) = (ml^sin^i)/ {nirg + m^d)^ (3) 

where mwd is the mass of the white dwarf, m^g is the mass of the red giant, and i is the orbital 
inchnation, to derive hmits on the masses and inchnation. Constraining the inchnation by 
the requirement that the observed jet is nearly in the plane of the sky required a very low 
mass white dwarf. The third star in the long period orbit was assumed to be an unseen 
late-type dwarf. The model presented in Paper I was based on a moderate inclination. As 
discussed in the present paper, this model is no longer tenable due to improved knowledge 
of the component masses and orbital inclination. 



5.3. Skopal et al. 1996 



Skopal et al.l (Il996a[ ) claimed that both long period and short period eclipses are present 
in optical photometr y and spectroscopy. They then proposed a revised three body model. 
Skopal et al.l (Il996bl ) felt that similarities between the 2.1 year period of CH Cyg and other 
symbiotic systems require that the symbiotic system be in the short period orbit. To produce 
the observable eclipses, the G-K dwarf in Paper I was replaced wit h a gian t in th e long period 
orbit. In recognition of the single lined nature of the spectrum, ISkopall (119971 ) argued that 
the spectrum of the giant in the long per iod orbit is blended with the M7 III primary, and 
one of the two giants is a bright giant. iTaranova &: Shenavriru (12004 ) suggested that the 
giants be nearly identical. 



This model has several problems. iMikkola fc Tanikawal (119981 ) found a triple system 
with a giant in the outer orbit to be unstable. There is no evidence in the spectra for the 
other bright giant. If two giants are present in the system, stellar evolution requires that 
they have nearly identical masses. In this case the inclination required by the mass function 
is relatively small and eclipses from the inner orbit will not occur. 



5.4. Mikkola & Tanikawa 1998 



Mikkola fc Tanikawal (119981 ) proposed that there is an inner close binary of total mass 
~4 Mq with a white dwarf of mass ^IMq in an outer orbit. The inner orbit is at a high 
inclination with respect to the white dwarf orbit. Activity on the white dwarf is driven by a 
Kozai resonance which causes large eccentricity variations in the inner binary. In the high- 
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eccentricity state, gas would be expelled from the red giant causing activity on the white 
dwarf. This model explains the long periods of inactivity in CH Cyg. 



The masses of the iMikkola &: Tanikawal (119981 ) model are not in agreement with the 
derived masses for the system, but the basic model is potentially viable. We will return to 
this point in discussing binary models for the short period variation. 



5.5. Two Star Models 



Munari et al.l (Il996l ) registered several objections to the triple star model for CH Cyg. 
These included orbital stability, proposed masses for the components, and modeling the 
photometric behavior. They suggested that t he short-period vari ations are pulsational, and 
hence, the system must be a binary. However, iMunari et al.l (119961 ) failed to note that a 2.1 yr 
pulsation exceeds the fundamental period. Thus, if the 2.1 yr period results from pulsation, 
it can be neither a normal pulsation mode of the star nor a beat between pulsation modes. 



Ezuka et al.l (119981 ) derived a minimum mass of 0.4 M© for the white dwarf, which 
exceeds the limit imposed in Paper I. They concluded that the value of the mass function 
determined in Paper I is implausible. While their minimum white dwarf mass is indeed 



inconsistent with the result for the 2.1 yr orbit. lEzuka et al.l (119981 ) failed to note that Paper 
I presented both short- and long-period mass functions, and the latter value is consistent 
with the 0.4 white dwarf mass limit. 



Schmidt et al.l (120061 ) briefly reviewed the controversy about whether the CH Cyg system 



is a triple or binary system. They also discussed whether the white dwarf is in the inner 
or outer binary. They noted the discovery that radial velocity variations of r nultiple-mode 
semiregular M-giants have periods longer than the fundamental radial mode (IHinkle et al. 
2OO2I ) . Indeed, CH Cyg has a strong similarity to these systems, which suggests that the 
2.1 yr period in CH Cyg could be a non-radial pulsation mode. If so, the long period 
velocity variation discussed in Paper I must correspond to the orbit containing the symbiotic 
components, and the CH Cyg system is a binary not a triple system. This is a viable model 
for the CH Cyg system, and we will examine this possibility more extensively below. 



6. THE CH CYG SYMBIOTIC SYSTEM 

Adopting sin^i = 1 and using the masses of 2.0 and 0.6 Mq that were estimated in 
§4.4, we compute from equation (3) a value for the mass function of 0.032 M©. This value 
is similar to our 15.6 yr mass function value of 0.051 Mq but quite different from our 2.1 
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yr mass function value of 0.0015 Mq. Increasing the whit e dwarf mass to 0.72 Mq, which 
imphes a progenitor mass of 3.3 Mq (IKahrai et aLll2008l ). produces a match to the 15.6 
yr mass function value of 0.051 Mq. This mass function value is then in agreement with 
multiple lines of evidence, reviewed above, that require the symbiotic components to be 
members of the long period system. For the above masses Kepler's third law produces a 
semimajor axis of 8.7 AU. Assuming that the 15.6 yr and 2.1 yr orbits are co-planar and 
that the red giant mass is 2.0 Mq, the 2.1 yr mass function requires a mass of ~0.2 Mq for 
a postulated secondary star in that orbit, a point that we will return to in discussing the 
nature of the short period variation. 

The masses of the dwarf and giant are, of course, related but not uniquely constrained 
by the 15.6 yr mass function. From the various limits set on the values for the masses, the 
red giant mass is in the range 1.5 < Mrg < 3.0. From the long period mass function the 
corresponding range for the white dwarf is 0.61 < M^^ < 0.92. Significantly, all the values 
for masses of the white dwarf derived from the long period mass function are in accord with 
evolutionary constraints. The mass range is sufficiently small to ex clude possible models fo r 
CH Cyg with either low or high mass white dwarfs (see for example iLuna fc Sokoloskil 120071 ). 



Iben &: Tutukovl (119961 ) commented in their extensive paper on the evolution of symbi- 



otic stars that symbiotic systems are a "nonuniform family of wide binaries with actively 
interacting components that differ in (a) the nature of their accreting components, (b) rea- 
sons for mass exchange, and (c) the physical mechanisms for their observed variability." This 
captures the difficulty in understanding CH Cyg. Although the system is phenomenologi- 
cally similar to many S-type symbiotics, in many respects it is a different type of symbiotic 
system. The orbital period of 15.6 yr is more than thre e times that of the next longest S-type 
symbiotic with a well defined orbit (IFekel et al.l 120071 ). In addition, none of the previously 
studied S-type symbiotics have a pulsating AGB star in the system. The D-type symbiotics 
R Aqr and o Cet are fundamental AGB pulsators (Miras), while CH Cyg is an overtone AGB 
pulsator. Consequently CH Cyg has a lower mass loss rate than the fundamental mode Mi- 
ras. However, in the CH Cyg system the lower mass flow to the secondary is compensated 
by the s maller semimajor axis . The much less studied symbiotic PU Vul is possibly a similar 
system ( iNussbaumer fc Vogellll996l ). 



The Mira symbiotics R Aqr and o Cet are similar to CH Cyg in some respects. For 
example, R Aqr has a jet, although its 40.9 year orbital perio d is nearly three times that 
of CH Cyg. The R Aqr system semimajor axis is 13-16 AU (IGromadzki fc Mikolajewska 
20081 ). o G et has a period m ore than ten times longer, ~500 yrs, with 70 AU component 
separation (jPrieur et al.l |2002| ) . o Cet exhibits flickering on the secondary with a time scale 
similar to that seen on CH Cyg (jWarnerlll972l ). Such flickering is not a common feature of 
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symbiotic systems. C H Cyg, R Aqr, and o Get have all been detected at X-ray wavelengths 
( iKarovska et al.ll2007l ). No doubt the comprehensive information about these three objects, 
including the X-ray flickering, is also a measure of the detail in which these nearby objects 
can be studied. X-ray flickering could well be present in other systems but is masked by 
other X-ray emission from the systems or by the larger distances to other systems. 

It is interesting to compare the near-IR colors of CH Cyg, both outside and during 
dust formation events, with the colors of symbiotic (D-type) Miras. The s pectr al type of 
CH Cyg is similar to that of many Miras near maximum light. IWhitelockl (119871 ) provides 
infrared photometry for a selection of symbiotic Miras. When quiescent, CH Cyg occupies 
the domain of normal Miras in a (J-K) - (K-L) relation. During a dust formation episode, 
the CH Cyg colors evolve alo i ig the line occupied by symbiotic Miras. Of the symbiotic 
systems studied by IWhitelockl (119871 ) CH Cyg is most similar to Hen 2-38. 



Whitelockl (119871 ) found it unlikely that symbiotic systems that currently contain Miras 



would have been recogn ized as symbiotic prior t o the onset of the Mira high mass-loss phase. 
She also noted, as did llben fc Tutukovl (Il996l ). that there is not an evolutionary relation 
between symbiotics. Symbiotics are binary systems of different separations, passing through 
an evolutionary phase that results in mass transfer. 

A question frequently raised about CH Cyg is the cause of the changes in the activity 
state of the symbiotic system. As discussed by iMikolajewski et al.l (Il990al ) CH Cyg was 
quiescent from the first observations in 1885 until 1963. The interval since 1963 has been 
dominated by periods of activity m ixed with quiescent periods. With a 15.6 year binary 
orbit and normal M giant mass loss, iMikolajewska et al.l (119881 ) found that the Bondi-Hoyle 
accretion onto a ~1 Mq white dwarf from the giant wind was lower by about an order of 
magnitude than the accretion required to drive the hot component luminosity. 

A solution is that widely separated, high mas s-loss rate symbiotic systen i s are under- 
going wind Roche-lobe overflow (see for instance iPodsiadlowski &: Mohamedl 120071 ). The 
stellar wind in these systems is directed by the Roche lobe and the efficiency of mass trans- 
fer to the secondary is much higher than might otherwise be expected. Models suggest mass 
transfer efficiencies approaching 100% instead of a few percent as previously believed. This 
ad dresses the active states but n ot the periods of inactivity. The latter have been attributed 
by ISokoloski fc KenyonI (l2003bl ) to the collapse of the inner accretion disk. We will discuss 
the activity of CH Cyg below in relation to the nature of the 2.1 yr variation. 



- 20 - 



THE 2.1 YR PERIOD 



In Paper I we assumed that both the 2.1 yr photometric and spectroscopic changes were 
due to orbital motion. With the recognition that the symbiotic system resides in the 15.6 yr 
orbit, we seek an understanding of the 2.1 yr periodic variation. 

Semiregular variables obey a period-luminosity relation ( Wood fc Sebol 1996h. From 
MACHO data a number of overtone series were identified in the LMC by lWood et al.l (119991 ). 
Since the publication of this seminal work, the relations have be en refined with data from 
other surveys, most recently those of OGLE (ISoszyhski et al.ll2007l ). Linear pulsation periods 
calculated by IWoodI (120071 ) can be compared with the LMC data to identify the modes of 
pulsation. Period-luminosity sequences are present for fundamental through fourth overtone 
pulsation. 

CH Cyg is well established as a semiregular variable with a period of ~100 days. Its 
absolute K ma gnitude is — 7 -5 ± 0.4 (§4.2.1). To compare CH Cyg with the LMC rela- 
tions found bv IWood et al.l (119991 ) and others, we take a distance modulus of 18.58 from 
Szewczyk et al. ( 20081 ) foi~the LMC. The pulsation mode(s) of CH Cyg then can be iden- 
tified from the appropriate period-luminosity relation. The CH Cyg 100 day period and 
magnitude are within 0.2 mag of the relation for first overton e pulsation (Fig. 4). This is 
well within the ~0.5 magnitude natural width of the relation (ISoszyhski et al.l 120071 ). 



In §4.1 the 2.1 yr period was identified as a "long secondary period" (LSP). LSP variables 
are so named because they all also have shorter period pulsation. LSP variables were found 
by lWood et al.l (119991 ) to obey a P-L relation (his sequence "D"). The CH Cyg 2.1 yr period 
is ~0.3 magnitude off the center of the D sequence (Fig. 4) but again well within the natural 
width of the relation. IWood et al.l (119991 ) pointed out that the range in the ratio of the long 
to the short period is a function of the long period and at per iods near 1000 days ranges 
from 5 to 15. The long to short period ratio for CH Cyg is ~8. iDerekas et al.l (120061 ) found 
that no normal stars fall between the fundamental and LSP sequences. Figure 4 illustrates 
that the 2.1 year CH Cyg period is clearly not fundamental pulsation. 



While approximate ly 25 - 30% of all pulsating AGB stars show LSP behavior (IWood et al 
19991 : iPercy et al.l 120041 ) , there is considerable disagre ement about t he ph y sical cause. Pos- 



sible o r igins of the LSPs h ave been investigate d by iHinkle et al.l (120021) , lOlivier fc Wood 



(|2003[ ). IWood et al.l (120041 ) . IPerekas et all (l2006f ). and ISoszvhskil (l2007n . among others. Six 
causes have been proposed: (1) radial pulsation, (2) non-radial pulsation, (3) a low mass 
companion, (4) rotating spheroidal shape for the star as the result of common-envelope evo- 
lution, (5) circumstellar dust clouds, and (6) star spots. Of this list, two causes of the LSPs 
remain viable after detailed investigations by the above authors: orbital motion from a close. 
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low-mass companion and non-radial pulsation. IWood et al.l (120041 ) concluded that the most 
likely explanation fo r the LSPs is a l ow de gree g mode in the outer radiative layers of the 
AGB star. However, ISoszyhski et al.l (120071 ) and others noted that the sequence D variables 
form an extension of the ellipsoidal variable sequence. This suggests that at least some LSPs 
result from extremely low-mass contact binaries. 



CH Cyg appears on the iPayne-GaposhkinI (Il954j ) and lHoukI (119631 ) lists of LSP variables. 
While these authors knew that the long to short period ratio was in the range 5 to 15, the 
discovery that AGB variab les have a period-luminosity relation did not occur for another 
25 years (iFeast et al.l Il989f). and the iden tification of the LSP period- luminosity relation 



followed 10 years later (IWood et al.l Il999l ). Thus, it is not surprising that the 2.1 year 
period in CH Cyg has not been extensively investigated as a possible LSP. However, it 
turns out that some of the same arguments that have been investigated as possible causes 
of LSP variability have beer i investigated for CH Cyg . Based on the specifics of the 770 
day photometric variability, iMikolajewski et al.l (119921 ) conclu ded that the 770 day perio d 
is not caused by va r iable extinction, as previously proposed by IMikolajewski et al.l (Il990al ). 
Mikolajewski et al.l (1l992l ) also dismissed rotation as a cause since the rotational velocity 
required is > 13 km s~^ and would produce noticeable line broad ening. In §4.3 the CH Cyg 



line width was shown to be slightly more than half this value. IMikolajewski et al.l (Il992l ) 
hypothesized that the 2.1 yr period might be the result of spot s related to cyclic c hanges in 
convective cells. This origin for LSP behavior was ruled out by IWood et al.l (120041 ). 



7.1. Binary Hypothesis 



Soszyhski et al.l (120041 ) and others show that the sequence of ellipsoidal variables (se- 



quence E) overlaps Wood's sequence D. The period-luminosity relation requires that binar y 
systems causing D sequence photometric variations be contact binaries (jPerekas et al.ll2006l ). 
Based on the shape of the light and velocity curves, it is apparent that the short period sys- 
tem in CH Cyg is not a s tandard contact binary . Vel ocity curves f or various contact binary 
systems can be found in ILu fc Rucinskil (119991 ) and IWoodI (120071 ). Contact binary radial 
velocity curv e s hav e a characteristically sinusoidal shape because their orbits are circular. 
Adams et al.l (120061 ) found the sequence E objects in their sample to be first ascent, giant 
branch objects with very small envelope rnasses, probably as the result of a common envelope 
mass ejection event. However, ISoszyhskil (120071 ) noted that if the Roche lobes are not full, 
it is possible to transfer mass driven by a stellar wind. In such a situation the eccentricity 
can be increased. Thus, it is of interest to discuss the binary option for CH Cyg. 



In our discussion of the symbiotic CH Cyg binary system, a number of arguments have 
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been used to assign mass, effective temperature, and radius to the M giant (§4). In addition, 
the inchnation of the 15.6 yr orbit is known from echpses to an uncertainty of less than 
a degree. We initially assume the simple case of coplanar 15.6 yr and 2.1 yr orbits. It 
is then possible to solve for the mass of a possible 2.1 yr companion. Assuming a 2 Mq 
primary, a putativ e 2.1 yr companion ^ yould have mass 0.2 Mq. Using the mass - spectral 
type calibration of iBaraffe fc Chabrierl (119961 ) the hypothesized low-mass companion is then 
an M3 dwarf with an effective temperature ~3200 K. This is so similar to the much more 
luminous M giant that such a dwarf would be undetectable spectroscopically. The Roche 
lobe radii of a 2 and 0.2 binary are 260 Rq and 96 Rq respectively in a circular 
orbit. However, the eccentricity e of the 2.1 yr orbit is 0.33. Thus, the Roche lobe will be 
constantly changing from (l+e)a for apastron to (l-e)a for periastron. The radius of the 
giant is ~280 R©, so such a system is Roche lobe filling, and so, at each periastron passage 
extreme mass loss would be expected. 

In Paper I, instead of a main sequence red dwarf, a white dwar/ companion of approxi- 
mately 0.2 Mq was hypothesized. We now argue that a 0.2 Mq white dwar/ companion can 
not exist. The presence of 2 Mq star in the CH Cyg system requires that any progenitor 
of a white dwarf was initially more r nassive than 2 Mq, . Such stars produce white dwarfs 
that are more massive than 0.6 Mq . (IKalirai et al.ll2008l ). However, we note that 0.2 Mq 
white dwarfs do exist; iLiebert et al.l (120041 ) reported the discovery of such an object. Low 
mass white dwarfs are believed to result from mass transfer and common envelope evolution 
with a degenerate companion. There seems no need to invoke this evolutionary path in the 
CH Cyg system. Even if such a white dwarf companion could exist, a contact system with a 
AGB star would be exceedingly active. In addition, the X-ray observations of CH Cyg result 
in mass limits which require a more massive white dwarf (§4.2.2). 



Soszyhski et al.l (120071 ) finds that the binary star explanation of the LSP requires the 
radius to semimajor axis ratio to be ~0.4. Kepler's third law then requires a mass of 6 
Mq for the dwarf. However, masses larger than 2 Mq are excluded if the stars in the sys- 
tem have equal ages. If the mass of the dwarf and the giant are equal the mass function 
requires an inclination of ^8 °. This in turn excludes the inclined orbit Kozai resonance 



model of iMikkola fc Tanikawal (119981 ). On the other hand, a coplanar result is also evi- 
dence that the triple mod el is unphysical because triple systems are seldom if ever coplanar 
jMuterspaugh et al.l[2008h . 



The 2.1 yr velocity variation (Fig. 3) is not sinusoidal. This provides one of the most 
powerful arguments against the 2.1 yr binary hypothesis. In fact, the eccentricity of 0.33 is 
moderately large. Polar views of the 15.6 yr and 2.1 yr orbits are shown in Fig. 5. ISoker 
(120001 ) provides formulae for computing the time scale for tidal synchronization. Because the 
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low mass companion does not have sufficient mass to spin up the en velope of the M gia nt, the 
synchronization time scale is the merger time scale, ~1000 years (IWood et al.ll2004l ). The 
time scale for orbital circularization is longer than the time scale for synchronized rotation. 
Merger in the synchronization time scale makes this argument mainly academic. However, 
the circularization time scale is also short, ~5000 year s, so all contact bin ary orbits of this 
type would be expected to be circular. As noted by IWood et al.l (120041 ). there are many 
known LSP systems, so eccentric orbits, even if possible for one object, could not possibly 
be the case for the entire class. 



Soszyhskil (120071 ) suggested that mass lost by the red giant in a binary system with a low 
mass star follows a spiral pattern producing the LSP variations. Such spiral mass loss could 
then affect the symbiotic binary system. However, as discussed in §4.5 infrared observations 
of CH Cyg require a dust temp erature which place s the d ust well outside the the 15.6 yr 
long period orbit. Furthermore, iMikolajewski et al.l (119921 ) find that the 770 day variation 
in CH Cyg, while color dependent, is inconsistent with reddening. 

These arguments provide convincing evidence that the 2.1 yr period variation of CH Cyg 
does not result from binary motion. For the entire LSP group iDerekas et al.l (120061 ) used 
arguments based on a multivariate test of the amplitude-luminosity relation to conclude that 
the LSP phenomenon is in general not a result of contact binaries. 



7.2. Pulsation 



While it has been shown theoretically that radial pulsation can not be a cause of the 
LSPs, this is still frequently hsted among the possible causes. In the case of CH Cyg several 
authors h ave suggested that the 770 day period is the result of radial pulsation of the M 
giant (see iMunari et al.lll996l ). We emphasize that a period of 770 days can not possibly be 
a radial pulsation. The frequencies of radial pulsations in AGB stars c an be calculated as 
a function of bolometric magnitude (IHughes fc Woodlll990l : IWoodI 120071 ). The fundamental 
period as a function of absolute K magnitude is plotted on Fig. 4. At the luminosity of 
CH Cyg, K = —7.5 mag, the radial fundamental period is ~250 days. If the origin of the 
LSP is pulsational, the pulsations must be non-radial in nature. 



Spectro scopic studie s of LS Ps have been carried out by lHinkle et al.l (120021 ). lOlivier fc Wood 



(2003 


), and 


Wood et al. 


(2004) 



First, the radial velocity variatio ns among LSP stars are ne arly identica l. The velocity phase 



curves for seven LSP stars from iHinkle et al.l (120021 ) and IWood et al.l (120041 ) are shown in 



Fig. 6. The position of these objects on the period-luminosity relation is shown in Fig. 4. 
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Although the original orbit for one of these stars, AF Cyg, had e = 0.08 ± 0.20, given the 
large uncertainty, we have recomputed the orbit with e = 0.3 as shown in Fig. 6 and found 
a very acceptable fit to the velocities, although the computed velocity curve is in antiphase 
with those of the other stars. The 2.1 yr velocity phase curve of CH Cyg matches those of 
the other LSP stars. The remarkable similarity of the "orbits" is inconsistent with both the 
expectation of random orbital orientation as well as the circularization constraints discussed 
above. The variables included in Fig. 6 are large visual amplitude LSP systems, so the 
similarity of velocity amplitudes could be a selection effect. 

Second, the LSP velocity variations totally dominate the shorter pulsation period in the 
LSP objects studied. For instance, in CH Cyg the overtone pulsations are not detectable in 
the velocities (§3) and result in only a low amplitude photometric variation (§4.1). Thus, 
LSP velocity changes are not easily differentiated from orbital motion for stars where the 
photometric variation is poorly known. This could be pervasive among suspected binary 
systems. A possible example is the velocity curv e of the X-ray, suspected neutron star 
symbiotic binary HD 154791 (IGalloway et al.ll2002l ). which matches the parameters of the 
LSP "orbits" of the stars in Fig. 6. 



As noted by IWood et al.l (120041 ). velocity curves with the shape shown in Fig. 6 can 
be reproduced by rotating prolate spheroids. The symmetry of such a spheroid requires a 
rotation period twice the period of the vel ocity curve . This period is in agreement with 8 
km s~^ measured line widths in CH Cyg. iKiss et al.l (120001 ) proposed a rotating pulsating 
ellipsoid model for the unusual light curves of two semiregular variables. A rotating prolate 
spheroid could result from a recent binary merger, but it is difficult to understand how 
an entire class of variables following a period-luminosity law could be the result of recent 
mergers. 

The observed properties of a rotating prolate spheroid can be generated by low order 
non-radial pulsation. In the case of non-radial pulsation a period-luminosity relation would 
exist. Non- radial pulsat i on ca n be prograde, i.e. in the direction of stellar rotation, or 



retrograde. lUnno et al.l (119891 ) stated that waves traveling in retrograde, as opposed to 



prograde, result in identical velocity curves but with the sign reversed. This could explain 
the antiphase velocity curve of AF Cyg in Fig. 6. 



Hatzed ( 1l996l ) provides models for some low 1 sectoral (1 = — m) modes. Unfortunately, 
the models are only for m = 2, 4, and 6. For these modes there are at minimum 4 sectors 
around the star and the velocity variations are nearly sinusoidal. The velocity amplitude 
decreases with increa sing m becaus e for higher m there are a larger number of velocity zones 
with cancelling sign (jHatzeslll996l ). The amplitu d e and asymmetry of CH Cyg suggest an 
1=1 mode. From the relations given by iHatzed (119961 ) the 1 = 1 mode pulsation velocity 
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amplitude is about half of the observed velocity amplitude. An interesting aspect of non- 
radial modes is that the vertical scale height is small. 

Additional iri s ight i nto the LSP problem comes from K to early M red giants. For these 



stars iHenry et al.l (I2OOOI ) reported the prese nce of short period, i.e. time scales of a few days 
to weeks, radial and non-radial pulsations. iHenry et al.l (120001 ) found non-radial pulsations 
restricted to the hot side of the coronal dividing line at about K2. For an early K giant 
the radial periods are typically days. The K giants have secondary periods with lengths of 
hundreds of days. There is considerable uncertainty about the origin of the LSPs for these 
stars, as for the M giant LSPs. As is the case for AGB LSPs, rotation with star spots, 
non-radial pulsations, and low m ass companions ha v e bee n discussed as explanations. In 
the case of the KO III star tt Her, iHatzes &: Cochran! (119991 ) found the pressure scale height 
for the photosphere to be a factor of 10 smaller than the radial wavelength of the 613 day 
pulsations. 



The hypothesis of non-radial modes in late-type giants has been impeded by the knowl 
edge that gravity modes are evanescent in convective regions. This mak es their postulated 



existe nce in both AGB and red giai its difficult to understand. Thus, IHatzes fc Cochran 
(119991 ) suggested r-mode oscillations (IWolfflll996l ). which arise in rotating, convective fluids, 
as an alternative. However, this alternative raises a different problem for AGB stars. In 
AGB stars the rotational period results in a pulsation period which is too long. For r-modes 
the pul sation frequ ency is equal to the rotation frequency for 1 = 1 and longer for higher 
modes (jWolfl ll996|) . In the case of CH Cyg th e rotation pe r iod is more than twice the 2.1 
year pulsation period. A solution proposed by IWood et al.l (120041 ) rehes on the point that 
LSP non-radial modes occur only in stars that also have radial pulsation. The proposal 
is that radial pulsation thickens the radiative layer above the convective layer allowing the 
propagation of g-modes. 



Wood et al.l (12004 ) discussed light and color variations in LSP stars. These can be com- 
pared with the variations in CH Cyg (§4.1) and predictions from stellar pulsation. The re- 
ported changes in spectral type are consistent with pulsation r elated changes in temperature. 
The types of variation reported by iMikolajewski et al.l (119921 ) are fully consistei it with vari 
ations for other LSP stars and with variations expected from stellar pulsation (IWood et al. 



A complicati on in underst a nding t he CH Cyg systera has been the report of short 
period eclipses by ISkopal et al.l (Il996al ). IWood et al.l (120041 ) concluded that for LSP stars 
the strength of Ha is a function of phase. They interpreted this as changing chromospheric 
activity. At some phases the chromosphere is absent, while at other phases it covers up to 
70% of the stellar surface. This same phenomenon has been reported for the K giant LSP 
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stars . For a Boo the strength of He I 10830 A is correlated with the LSP ( Hatzes fc Cochran 
19931 ). W e suggest tha t the U V cha nges in CH Cyg co rrelated with the short period phase 
found by ISkopal et al.l (Il996al ) and lEyres et al.l (120021 ) are pulsat ion driv e n cha nges in the 
chromosphere. From the dates of the deep UV minima given in ISkopall (119951 ) the 2.1 yr 
orbital phase (from the eccentric orbit) is 0.14, which is conjunction in an orbital model or 
a time of null velocity movement in a pulsational model. The absence of flickering during 
at least some of these events could be un related changes of structure in the white dwarf 
accretion disk (ISokoloski fc Kenyonll2003al ). 



In Paper I the change of line shape with short period phase was discussed. The spectral 
lines have the largest FWHM at = 0.4-0.6. However, the lines are mo re narrow and 
symmetric at phase 0.4 becoming more asymmetric at phase 0.8. For K giants iHatzed (119961 ) 
has noted the change in bisector shape with phase and mode in non-linear pulsations. The 
resolution and signal-to-noise ratio of the current CH Cyg spectra does not allow detailed 
comparison but indicates that further study would be of interest. 

A long standing question for CH Cyg is why this star was an M giant spectral type 
standard for the first half of the 20th century and then became an active symbiotic star. At 
the heart of this question is the the nature of the accretion onto the white dwarf secondary. 
A possibly naive expectation is that the LSP variation is a controlling factor in the M giant 
mass loss. We speculate that in the case of a non-radial pulsation origin for LSP, the non- 
radial modes are stable but variations in the dominant mode can occur. If the dominant 
mode controls other stellar properties, for instance the mass loss rate or the directionality 
of the mass loss, the mass transfer in a symbiotic system would change. Observational tests 
employing samples of LSP variables should be possible. 



8. CONCLUSIONS 

The near-infrared radial velocities of CH Cyg conclusively show that the M giant in the 
system exhibits two different velocity variations, a 15.6 year "long" period and a 2.1 year 
"short" period. We have reviewed the literature on the basic parameters for the CH Cyg 
giant and have estimated its mass and radius. Evolutionary arguments require an M giant 
mass near 2 Mq. X-ray observations require a white dwarf companion more massive than 
0.44 Mq. Various lines of evidence but most compellingly observations of eclipses demand 
that the inclination of the 15.6 year orbit is nearly edge on. The long period mass function 
is then matched by the assumption of a 0.7 Mq white dwarf. 

The CH Cyg symbiotic system is an unusual one. The separation between the giant 



-27- 



and white dwarf is about four times larger than that for typical S-type systems. CH Cyg 
is a first overtone pulsating AGB star, which will drive symbiotic mass transfer at a larger 
radius. While CH Cyg is classified as an S-type system, there are a number of similarities 
with D-type Mira systems. Changing activity in the CH Cyg symbiotic system is possibly 
related to variations in the short-period activity. 

The observed 2.1 yr velocity variation is indistinguishable from that seen in stars with 
long secondar y periods (LSP). LS P objects occupy a track on the period-luminosity diagram, 
identified by IWood et al.l (119991 ) as sequence D. The cause of the LSP variations, while 
discussed exhaust ively, remaiiis unc ertain. LSP is the only theoretically unexplained type of 
stellar variability (IWood et al.ll2004l ). In the case of CH Cyg there seem to be only two viable 
options to produce the 2.1 yr period. The star could be a rotating prolate spheroid, perhaps 
as the result of a common envelope phase with a former low-mass companion. However, with 
this model it is hard to understand why there are so many LSP stars. Much more likely, CH 
Cyg is undergoing low order non-radial pulsations that mimic the observed properties of a 
rotating prolate shape. If these are g-modes, the outer stellar structure of pulsating M giants 
does not follow standard models because an outer radiative layer is required to propagate 
the g-modes. While this structure change is not required for r-modes, the periods of r-modes 
are excessively long. 
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Table 1 — Continued 
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Table 2. ORBITAL ELEMENTS AND RELATED PARAMETERS OF CH CYG 





Short-Period 


Long-Period 


Paranieter 


"Orbit'^ 


Orl)it 


P (days) 


750.1 ± 1.3 


5689.2 ± 47.0 


P (years) 


2.0537 ± 0.0036 


15.58 ± 0.13 


T (HJD) 


2,447,293.5 ± 12.9 


2,445,681 ± 192 


7 (km s"-*") 




-59.91 ± 0.09 


K (km s-i) 


2.87 ± 0.13 


4.45 ± 0.12 


e 


0.330 ± 0.041 


0.122 ± 0.024 


uj (deg) 


229.5 ± 7.7 


216.9 ± 12.7 


a sin i (10^ km) 


27.90 ± 12.30 


345.69 ± 9.09 


/(m) (Me) 


0.0015 ± 0.0002 


0.051 ± 0.004 
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Table 3. DERIVED CH CYG COMPONENT PARAMETERS 

Parameter Value Reference 

Red Giant: 

M^i -4.5 §4.2.1 

L 5000 Lq §4.2.1 

T 3100 K §4.2.1 

R 280 Rq §4.2.1 

M 2 Mq §4.4 

White Dwarf: 

L 0.25 Lq §4.2.2 

M >0.56 Mo §4.4 

Long Period Orbit: 

P 15.6 yr §3 

a 8.5 AU §4.6 

i 84° §4.6 

Circumstellar Shell: 

■dinner 22 §4.5 
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Fig. 1. — Our 106 KPNO radial velocities, derived from lines in the near infrared, of the 
CH Cyg M giant component plotted versus heliocentric Julian day. The data span nearly 30 
years. A short period, 2.1 year, and a long period, 15.6 year, variation are apparent. 
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Fig. 2. — The computed velocity curve of the 15.6 yr long period orbit compared with the 
KPNO radial velocities. Each plotted velocity consists of the total observed velocity minus 
its calculated short period velocity. Zero phase is a time of periastron. 
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Fig. 3. — The computed velocity curve of the 2.1 yr velocity variation interpreted as an 
orbit, compared with the KPNO radial velocities. Each plotted velocity consists of the total 
observed velocity minus its calculated long period velocity. Zero phase is a time of periastron. 
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Fig. 4. — Period-absolute K magnitude plot of field LSP stars. Each LSP star has two 
periods, a normal pulsa tion mode a n d the LSP. Filled circle is CH Cyg. Open squares 
are six LSP sta rs from iHinkle et al.l (120021 ) and IWood et al.l (120041 ) that have Hipparcos 
parallaxes. The Ivan LeeuwenI (120071 ) revised Hipparcos parallax values are used. The lines 
Pq, Pi, and P2 show respectively l inear fundamental, first overtone, and second overtone 
pulsation periods from lWood et al.l (1 19991) for oxveen-ric h opacities. The hne D is the middle 



range of the Wood's D relation (see IWood et al.lll999l ). For field stars the uncertainty in 



the distance, indicated by the magnitude error bars, combined with the intrinsic range in 
magnitude blurs the observed pulsation relations, although most of these stars are clearly 
first overtone pulsators. 
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Fig. 5. — Polar view, to scale, of the computed orbits of the CH Cyg system. On the left, 
the inner and outer ellipses are the 15.6 yr orbits with a red giant of mass 2 Mq and a white 
dwarf of mass 0.72 Mq (see text). The shaded circle and the dot represent the giant and 
white dwarf, respectively, at periastron. The shaded circle is the scaled size (R=280 Rq) of 
the red giant which is assumed spherical. The plus is the center of mass of the orbit. The 
dashed line box is shown enlarged on the right. The inner and outer ellipses are the 2.1 yr 
orbits of the red giant and the low mass companion shown at periastron. The plus sign is 
the center of mass of the 2.1 yr orbit. The 2.1 yr and 15.6 yr orbits are shown coplanar (see 
text). The companion object in the 2.1 yr orbit can not be a white dwarf. However, we 
argue (see text) that the 2.1 year orbit is not physical. 
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Fig. 6. — A phase plot of the calculated radial velocity "orbits" of long secondary periods 
for eight stars plus the 2.1 yr per iod of CH Cyg. Six stars, RS CrB, AF Cyg, X Her, g Her, 
V574 Oph a r id BI Peg are from iHinkle et al.l (120021 ). and two, Z Eri and S Lep, are from 
Wood et al.l (120041 ). The center of mass velocity of each "orbit" has been subtracted from 
the observations to produce velocity curves with the same zero point for ease of comparison. 
The "orbit" of CH Cyg is a solid line, the other stars are dashed lines except for AF Cyg, 
which is a dotted line to emphasize that its "orbit" is in antiphase with the others. 



